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Abstract. We present BVI photometry and long-slit Ha rotation curve data obtained with ESO VUr/F0RS2 for six low 
surface brightness galaxies with extremely blue colours and very faint central regions. We find no evidence for a steep central 
density cusp of the type predicted by many N-body simulations of cold dark matter (CDM) halos. Our observations are instead 
consistent with dark matter halos characterized by cores of roughly constant density, in agreement with previous investigations. 
While unremarkable in terms of the central density slope, these galaxies appear very challenging for existing CDM halo models 
in terms of average central halo density, as measured by the Av/2 parameter. Since most of our target galaxies are bulgeless 
disks, our observations also disfavour a recently suggested mechanism for lowering the central mass concentration of the halo 
by means of a fast collapse phase, as this scenario predicts that the original CDM profile should still be detectable in bulgeless 
galaxies. Other potential ways of reconciling the CDM predictions with these observations are discussed. 
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1. Introduction 

The cold dark matter (CDM) scenario, in which the dark mat- 
, ter particles are assumed to be non-relativistic at the time of 
■ decoupling, and to interact predominantly through gravity, has 
been very successful in explaining the formation of large-scale 
structures in the Universe (e.g. Primack 2003 1. On the scales of 
galaxies, the CDM predictions of halo shapes, substructure and 
density profiles have however not yet been confirmed in any 
convincing way. 

Dwarf galaxies and low surface brightness galaxies 
(LSBGs) are believed to be more or less completely dominated 
by dark matter, thereby making them among the best probes of 
the density profiles of dark matter halos on galactic scales. A 
lot of recent research in this field has revolved around apparent 
discrepancies between the CDM halo density profiles predicted 
by N-body simulations (e.g. Navarro et al. 1996 1997; here- 
after NFW) and the dark halo density profiles inferred from 
observations (see e.g. de Blok et al. 2 001,. de Blok & Bosma 
[2002 1. Both the observed slope of the innermost density profile 
(the core/cusp problem) and the overall shape of LSBG rota- 
tion curves are in conflict with the N-body results, indicating a 
potentially serious problem for the CDM scenario. 
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To remedy this situation, numerous solutions have been 
proposed. Some of the discrepancies may be reduced by as- 
suming that dark matter is not cold, but rather self-interacting 
(Spergel & Steinhardt .2000j . warm (Bode et al. i2001i) . anni- 
hilating (Kaplinghat et al. l2000L decaying (Cen l2001> . or fliat 
the dark energy of the Universe has a phantom-like equation of 
state (Kuhlen et al. 2005 1. Another option may be to drop the 
notion of dark matter altogether and instead modify the laws of 
gravity (e.g. McGaugh & de Blok 1998b i. Most observational 
studies so far have however assumed the dark matter halos to 
be spherical, whereas CDM in fact predicts triaxial halos. It 
has been suggested that a realistic treatment of disk dynamics 
inside these more complicated potentials may possibly remove 
the discrepancy (Hayashi et al. l2004bV Yet another possibiUty 
is that the dark matter domination of some target galaxies have 
been overestimated (Graham 2002, Fuchs l2003l and that some 
complicated baryonic process may have significantly affected 
the density profile. 

Despite a large number of failures to confirm the NFW pre- 
dictions for the CDM halo profile on galactic scales, there are at 
least two objects for which the data have been reported to be of 
sufficient quality to allow well-constrained fits, yet NFW pro- 
files still appear reasonably consistent with the observations: 
NGC 5963 (Simon et al. 2005 1 and DDO 9 (de Blok 2005). Is 
it then possible that some galaxies have NFW halos whereas 
others do not? Could it be that the apparent discrepancy be- 
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Fig. 1. Isophotal B-band images of the target galaxies. The field 
size is 2' x2', except for ESO 146-14, for which it is 4' x4'. The 
outermost isophotes plotted correspond to //^ = 24.5 + 0.2 mag 
arcsec"^, and the flux steps between the brighter contours are 
constant on a linear scale (but different from galaxy to galaxy, 
to clearly display the structures present). North is up and east is 
to the left. The horizontal line visible in several of the images 
is caused by the gap between the two F0RS2 CCDs . 



tween predicted and observed halos is due to some bias related 
to how the target galaxies have been selected? 

Here, we investigate the central kinematics of a class of 
LSBGs not previously targeted by similar investigations. Our 
targets are bulgeless disks with extremely blue colours and very 
low surface brightness centres; properties which could possibly 
imply more dark-matter dominated central regions and allevi- 
ate potential problems of baryon-domination. In Sect. 2 we de- 
scribe the observations and the selection criteria of our galaxy 
sample. In Sect. 3 the rotation curves are presented and in Sect. 
4 the luminosity profiles. Sects. 5 and 6 compares the observed 
average central density and the inner slope of the density profile 
to the CDM predictions. Sect. 7 discusses the implications of 
our results and a number of potential problems with our analy- 
sis. In Sect. 8, our findings are summarized. 



Table 1. List of objects, integration times and observation 
dates. Photometric observations are labeled B,V or / depend- 
ing on filter. Spectroscopic long-slit observations are labeled = 
and II for slits aligned along the major and minor axis of the 
galaxy, respectively. 



ESO id. 


Integration time (s) 


Date (yymmdd) 


031-13 


B(852) 


020801, 020813 




V(372) 


020801, 020813 




7(1040) 


020801, 020813 




=(6324) 


020913 




11(1563) 


020831 


146-14 


B(120) 


020706 




K(68) 


020706 




/(1 20) 


020706 




=(2664) 


020718 


462-32 


B(320) 


020506 




K(136) 


020506 




7(372) 


020506 




=(6324) 


020518 




11(1560) 


020518 


532-32 


B(270) 


020706 




V(120) 


020706 




/(304) 


020706 




=(3162) 


020714 




11(1560) 


020706 


546-34 


B(744) 


020711, 020801 




V(340) 


020711, 020801 




7(880) 


020711, 020801 




=(6324) 


020911 




11(2080) 


020815 


548-09 


B(542) 


020801 




V(256) 


020801 




7(744) 


020801 




= (6324) 


021005 




11(2080) 


021005 



2. Observations, data reduction and sample 
properties 

2.1. Selection criteria 

The six disk-like LSBGs used in this investigation were se- 
lected from the ESO-Uppsala catalogue (Lauberts & Valentijn 
I1989> with the criteria that the targets should have: 

1) An average B - R < 0.5 mag inside the region with B- 
band surface brightness 20.5 < fig (mag arcsec"^) < 26; 

2) A B-band surface brightness inside the central 5" of 
fj.B,Q S 22.3 mag arcsec"^; 

3) A high inclination, but be separated from edge-on pro- 
jection by at least 5° (under the approximation of an infinitely 
thin disk). 

The first and second selection criteria are similar to those 
used in our previous studies of extremely blue LSBGs (e.g. 



Ronnback & BergvaU [T994I Bergvall et al. lT999l Zackrisson 
et al. 2005i, which indicate that objects of this type are metal- 
poor and almost extinction-free. The galaxies selected are ESO 
031-13, ESO 146-14, ESO 462-32, ESO 532-32, ESO 546-34 
and ESO 548-09. Out of these, ESO 146-14, 546-34 and 548- 
09 have featured in our previous papers on blue LSBGs. As far 
as we can tell, LSBGs of this type have not been targeted by 
previous dark halo investigations using optical rotation curve 
data (e.g. de Blok et al. 2001 , de Blok & Bosma 2002, Swaters 
et al. 2003 1. The selection criteria are further discussed in Sect. 
2.3. 

2.2. Observations and reductions 

All observations were carried out in service mode with the 
F0RS2 instrument at die VLT-UT4 Yepun 8.2m ESO telescope 
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Table 2. Target galaxy data. Here, and Mg represent the apparent and absolute integrated fi-band magnitudes inside the fig o 
- 26.5 mag arcsec"^ isophote. B - V and V - I denote the integrated colours inside the same radius, whereas <B - R> (taken 
from the ESO-Uppsala catalogue; Lauberts & Valentiin 19891 is the average colour inside the 20.5 < fig (mag arcsec"^) < 26 
region, /i/ represents the /-band scale length derived from the outer part of the surface brightness profile. Both the true B-band 
central surface brightness, figxh and the B-band surface brightness of an exponential profile extrapolated to the centre, have 
been corrected for inclination, /, assuming an infinitely thin disk. Vsys.hei denotes the systemic, heliocentric velocities and D the 
estimated distances. All magnitudes have been corrected for galactic extinction (Schlegel et al. 1998 1. 



ESO id. 


l'sys,hel 


D 


mjj 


Ms 






B-V 


<B-R> 


V -I 


h, 






(km s-') 


(Mpc) 


(mag) 


(mag) 


(mag 


(mag 


(mag) 


(mag) 


(mag) 


(kpc) 


D 












arcsec"^) 


arcsec"^) 












031-13 


6401 


82.7 


16.7 


-17.9 


24.3 


23.2 


0.49 


0.28 


0.95 


3.5 


83 


146-14 


1686 


21.0 


14.9 


-16.7 


23.5 


24.0 


0.21 


0.25 


0.70 


2.8 


82 


462-32 


2848 


38.7 


16.7 


-16.3 


23.8 


23.2 


0.45 


0.38 


1.31 


1.5 


83 


532-32 


2674 


35.8 


16.3 


-16.4 


24.4 


23.4 


0.42 


0.48 


0.60 


2.4 


79 


546-34 


1582 


19.5 


15.5 


-15.9 


24.2 


22.8 


0.21 


0.41 


0.54 


1.3 


76 


548-09 


1832 


22.7 


17.3 


-14.4 


24.7 


24.4 


0.45 


0.38 


0.80 


1.5 


80 



in 2002. Imaging in BVI was carried out at a seeing of < 1.2", 
while long-slit spectroscopy along both the major and minor 
axis of the target galaxies was carried out at a seeing of < 0.8". 
The spectroscopy was performed using the 1.0" slit and the 
1200R grism, with a wavelength range of 5750-7310 A and a 
spectral resolution of 35 km s"' per pixel at the central wave- 
length. Contour plots of the B-band images are displayed in 

Fig.m 

Observing dates and exposure times are summarized in 
Table. Due to technical problems, the minor-axis obser- 
vations for ESO 146-14 were unfortunately not properly ex- 
ecuted. The pipeline-reduced images were recalibrated us- 
ing standard stars, stacked and sky subtracted using the 
ESO-MIDAS package. Spectra were reduced and wavelength- 
caUbrated from arcspectra, using the same software. After 
stacking, all spectra were recalibrated to spatial bins of 0.75", 
except for the major-axis data of ESO 532-32, where the low 
signal-to-noise required the use of bins 1.25" wide. 

Galaxy centres and inclinations were determined from the 
outer isophotes of the images, the latter under the assumption of 
an infinitely thin disk. Surface brightness profiles were derived 
from in-house software developed exclusively for this purpose. 
Distances to all galaxies were estimated from the systemic, he- 
liocentric redshifts after converting to the centroid of the Local 
Group using lAU specifications and correcting for Virgocentric 
infall using the model by Schechter il980> . A distance to the 
Virgo cluster of 17 Mpc, corresponding to a Hubble constant 
of Ho = 75 km s"' Mpc"', was assumed. All magnitudes were 
corrected for Galactic extinction using the B-band extinction 
maps by Schlegel et al. ( 1 998 ), and the standard (R v = 3.1) ex- 
tinction law (Cardelli et al. 1989 1, as implemented in the NED. 
Due to the different extinction corrections used here, the data 
for ESO 146-14 and ESO 546-34 differ slightly from those pre- 
sented in Zackrisson et al. J2005t . 

Table 121 summarizes the BVI magnitudes, heliocentric sys- 
temic velocities, distances, disk scale lengths, inclinations and 
central surface brightness levels of our objects. Here, fisfi refers 



to the true B-band central surface brightness (integrated over 
the central 1.2", which corresponds to the maximum seeing 
disk), whereas fj^ q refers to the surface brightness of a fitted ex- 
ponential disk extrapolated to the centre. The disk scale lengths 
are derived from the outer parts of the /-band surface brightness 
profile, as discussed in Sect. 4. 

2.3. General sample properties 

For measuring the dark halo density profile, it would appear ad- 
vantageous to target the most dark-matter dominated galaxies, 
thereby minimizing the complicated effects associated with lu- 
minous baryons. Since the relation between galaxies and their 
dark halos is far from well-understood, this exercise is however 
not trivial. 

The galaxies used in this investigation have been selected 
on the basis of faint central regions, high inclinations and very 
blue colours. What bearing may these selection criteria have 
on the dark matter properties of these objects and the accuracy 
with which they may be derived? 

In McGaugh & de Blok ( 1998a) , a strong correlation was 
demonstrated between the dynamical mass-to-light ratio and 
the central disk surface brightness. This relation would im- 
ply that selecting the LSBGs with the faintest central regions 
should be advantageous for testing CDM models, unless the 
baryonic mass-to-light ratio is also unusually large for these ob- 
jects. To target the most extreme low surface brightness galax- 
ies, we have selected objects with a surface brightness inside 
the central 5" of /io,B > 22.3 mag arcsec"^, in combination 
with a high inclination. After correcting for inclination, four 
out of our six objects turn out to have yUo.B > 24 mag arcsec"^, 
making them among the most low-surface brightness objects 
used so far for measuring the density profiles of dark halos. 

Five out of our six targets (the exception being ESO 146- 
14) are underluminous in the central region compared to what 
would be expected from an exponential disk fitted to the outer 
isophotes. Hence, they appear to be completely bulgeless ob- 
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jects, which should make them ideal for testing the scenario 
advanced by Mo & Mao J20U4> to explain the discrepancy be- 
tween observations and the halo profiles predicted by CDM. 
Mo & Mao suggest that feedback associated with a fast col- 
lapse phase of baryons during the galaxy formation process 
could lower the high central densities predicted by CDM and si- 
multaneously explain the presence of bulges. If this scenario is 
correct, the original CDM halo profiles should therefore still be 
observable in bulgeless disk galaxies dominated by dark mat- 
ter, such as the ones presented here. 

Selecting galaxies with high inclinations comes with ad- 
vantages as well as disadvantages. By targeting disk galaxies 
close to edge-on, it is possible to identify the most extreme 
LSBGs in surveys limited by surface brightness, and to mini- 
mize the inclination correction to line-of-sight velocities mea- 
sured by long-slit spectroscopy, as well as the uncertainties in 
the position angle. High-inclination galaxies are however also 
sensitive to extinction and projection effects, especially in the 
central regions. The exact inclination limit at which these prob- 
lems become severe is unfortunately somewhat uncertain, as 
further discussed in Sect. 7.1. 

Since blue colours imply a lower stellar mass-to-light ratio 
M/L* (e.g. Bell & de Jong 2001 1, one may naively expect that 
blue LSBGs should be more dark-matter dominated than their 
red counterparts for a given central luminosity (surface bright- 
ness). This conjecture does however assume that blue and red 
galaxies are located inside identical dark matter halos, which 
is not necessarily the case. Graham ( 2002 1 examined the ratio 
of stellar disk to dynamical mass and did indeed find it to de- 
crease for bluer galaxies. Zavala et al. (■2003i) on the other hand 
estimated the ratio of baryonic to dynamical mass and found 
a possible trend in the opposite direction. From an evolution- 
ary point of view, the Zavala et al. result can be understood in 
a scenario in which the dark halos that formed early are more 
concentrated than those that formed late (Navarro et al. 1997 
Bullock et al. 2001 1, and in which blue galaxies are younger 
than red ones in an absolute sense. Although the last condition 
may seem reasonable, it is difficult to prove for LSBGs at the 
present time (Zackrisson et al. 2005 1, since the star formation 
histories may also differ between red and blue stellar popula- 
tions. 

Although the relation between colour and dark matter prop- 
erties is obscure at best, targeting very blue LSBGs does have 
other advantages. Since these objects are seen to be low- 
metallicity objects with little internal extinction (e.g. Ronnback 
& Bergvall 1994 Bergvall et al. 1999 1, chemical evolution is 
less of a concern when estimating the Mj of their stellar pop- 
ulations (Zackrisson et al. 2005 ). When retrieving the rotation 
curve from high-inclination disks, dust effects are also likely 
to be minimized. We may furthermore be confident that the 
central surface brightness levels of our targets have not been 
significantly underestimated due to dust effects. 

3. Rotation curves 

Rotation curves have been derived along both minor and major 
axes by means of fitting a gaussian to the line profile of the Ha 
(6563 A) emission Une. For edge-on disks, this fitting proce- 
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Fig. 2. Position-velocity diagram for the Ha data of ESO 03 1- 
13. The contours trace the Ha flux on a logarithmic scale, 
whereas the dots indicate the rotational velocities resulting 
from a Gaussian fit to the line profile in each radial bin. The 
dashed vertical line indicates the adopted centre of the galaxy 
and the dashed horizontal line the systemic velocity. The line 
profiles do not display any obvious systematic tails towards the 
systemic velocity. 

dure may underestimate the true velocity because of projection 
effects, since a substantial amount of the light can originate 
at greater radii in the disk, where the line-of-sight velocities 
are lower In these cases, more sophisticated fitting techniques 
should be employed to recover the true rotation velocity (e.g. 
Kregel & van der Kruit 2004 Gentile et al. 2004 1. Although the 
inclinations of our targets are high and based on the approxima- 
tion of an infinitely thin disc (therefore probably slightly under- 
estimated), the emission-line profiles of these galaxies do not 
show any obvious systematic tail towards the systemic velocity, 
as would be expected if projection effects were significant. As 
an example, the Ha position-velocity diagram for ESO 031-13 
is displayed in Fig. |2l The lack of systematic asymmetries in 
the line profiles leads us to assume that the disks are very thin 
and that kinematic projection effects are small. The possibility 
that a slight systematic skewness could be masked by the lim- 
ited spectral resolution of these observations, in combination 
with a patchy distribution of emission-line regions in the disk, 
can however not be completely ruled out. See Sect. 7.1 for a 
more detailed discussion about potential projection effects. 

The approximate centre of each galaxy was estimated from 
the outer isophotes of the continuum radiation. Because of the 
flat light profile encountered in many cases (see Sect. 4), the 
exact centre was then determined by maximizing the symme- 
try between the two sides of the rotation curve inside a limited 
spatial region. In order to derive mass profiles of the galaxies, 
the major-axis rotation curve was folded and spatially averaged 
over bins large enough to produce a reasonably smooth appear- 
ance. No artificial smoothing (of the type used in e.g. de Blok 
et al. 12001' and de Blok & Bosma 2002 1 has however been em- 
ployed. The error bars on the spatially averaged rotation curves 
were simply taken to be the uncertainties in the mean veloc- 
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Fig. 3. The major-axis rotation curves of the target galaxies. Black and white filled circles represent data from approaching and 
receding sides, respectively. The exception is ESO 546-34, for which the amplitude of the rotation curve is too low and the 
velocity field too disturbed to accurately determine this. For all other objects, symmetrised and spatially averaged rotation curves 
are plotted as black lines. For each galaxy, a dotted, vertical line indicates the radius at which the surface brightness profile starts 
to deviate from the exponential disk determined from the outer isophotes. 



ity of each bin. A minimum uncertainty of 5 km/s was finally 
imposed, to avoid unrealistically small errors associated with 
regions of the spectra with very high S/N. 

Both the raw and the symmetrised, spatially averaged 
major-axis rotation curves are displayed in Fig. |3] The corre- 
sponding minor-axis velocity profiles are shown in Fig.0] The 
latter have been folded along the spatial axis only (not in ve- 
locity). Linear distances along the radial axis of the minor-axis 
velocity profiles have been corrected for inclination assuming 
an infinitely thin disk. Since these corrections are very large at 
high inclinations, the linear scale is admittedly highly uncer- 
tain, and in many cases in poor agreement with the scale of the 
major axis rotation curves. 

A component of minor-axis rotation, as that found by 
Hayashi et al. f2004b') in their simulation of a disk rotating in a 
triaxial CDM halo, would in Fig.0]turn up as velocities of op- 
posite signs for the two sides. In no case is there any clear-cut 
evidence for this efifect among our target objects. 

In this paper, we make no detailed attempts to separate the 
dififerent contributions (e.g. stellar disk, HI disk, dark halo) to 
the observed rotation curves, as HI rotation curves are currently 
lacking for these galaxies. While interesting constraints on the 
spatial distributions and relative masses of the different com- 
ponents may sometimes be inferred even without HI data (e.g. 
Borriello & Salucci .2001j . the disc-halo degeneracies (e.g. van 
Albada et al. I1985> in the region covered by optical rotation 



curves are often severe. An attempt to perform a mass decom- 
position of these galaxies will however be presented in a future 
paper (Mattsson et al., in preparation). 

3.1. Notes on individual galaxies 

Because of the high resolution of Ha rotation curves and the 
typically patchy distribution of star-forming regions in LSBGs, 
irregularities of varying severity are present in all rotation 
curves presented here (although probably not more so than 
in previous investigations - see e.g. de Blok & Bosma 1^0021 
Swaters et al. '20031 In a couple of cases, the rotation curves 
are so disturbed that there is little point in trying to uncover the 
underlying mass distribution from the observed kinematics. 

ESO 031-13: The rotation curve is regular, except for at 
small wiggles at 3-6 kpc from the centre. The minor axis veloc- 
ity profile shows no signs of significant non-circular motions. 

ESO 146-14: The rotation curve is full of wiggles and 
asymmetries. Indications of this were already present in 
the lower-resolution rotation curve of Bergvall & Ronnback 
(1995). Because of these small-scale irregularitites, a sym- 
metrised, average rotation curve (admittedly highly uncertain) 
is constructed from very large radial bins. No minor-axis data 
is available (see Sect. 2.2). 

ESO 462-32: The rotation curve shows signs of non- 
circular motions in the centre, but is otherwise regular. The 
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ESO 031-13 



ESO 462-32 



ESO 532-32 




Fig. 4. The minor-axis velocity profiles of the target galaxies. Black and white filled circles represent data from opposing sides of 
the centre. The rotation curves have been folded spatially, but not in velocity. The radial axis has been corrected for inclination. 



minor axis kinematics are unusual, with redshifted features on 
both sides of the centre. 

ESO 532-32: While the major-axis rotation curve looks 
reasonably regular, the minor-axis velocity field does indicate 
the presence of non-circular motions at large distances from the 
centre. 

ESO 546-34: The major-axis velocity profile looks ex- 
tremely strange with very little net rotation. The minor-axis 
data also indicates strong non-circular velocities. Although the 
outer isophotes of this galaxy appear fairly regular, there is an 
indication of a warp-like feature in the centre (see Fig.Q, rem- 
iniscent of two disks partly overlapping each other It therefore 
seems likely that this is really two galaxies in a late stage of 
merging. Since no reliable mass profile can be derived from 
these data, no symmetrised and rebinned rotation curve is pre- 
sented. 

ESO 548-09 The rotation curve displays significant wig- 
gles (seen only on one side of the galaxy at a time) at around 
1 .3 and 2.3 kpc from the centre. The minor-axis velocity profile 
is consistent with no significant non-circular motions. 

4. Luminosity profiles 

The /-band surface brightness profiles, depicted in Fig. |5] are 
derived by integrating the light over elliptical strips of constant 
orientation determined from the inclinations and position an- 
gles of the outer isophotes. 

All galaxies except ESO 146-14 display a central light de- 
pression relative to the expectations from an exponential disk 



profile fitted to the outer parts of the galaxy. The B-V and V-I 
colour profiles typically do not show any dramatic gradients in 
the region when the depression starts to develop, indicating that 
this feature is unlikely to be an effect of dust reddening. The 
only exception is ESO 031-13, where the /-band profile is sig- 
nificantly steeper in the centre than B or V. This V-I gradient 
is however visible outside the break in the surface brightness 
profile as well. Since there is no corresponding slope in B -V, 
this is more likely to be due to a gradient in stellar population 
properties rather than to dust. In general, the colours vary very 
little across the face of the disks. 

The reason for the central depression in the surface bright- 
ness profile of LSBGs is not well-understood. These profiles, 
which are seen in many blue LSBGs (e.g. Ronnback & Bergvall 
1994 Bell et al. 2000 1 differ from the Freeman Type II pro- 
files (Freeman 1970J seen in high surface brightness galaxies 
(e.g. MacArthur et al. |2D03 1 in that the decrease in surface 
brightness persists all the way to the centre. Since no signif- 
icant central colour gradients are seen in the optical/near-IR 
data (Bergvall et al. 1999 (, these features cannot in general be 
attributed to dust eff'ects. If this central light depression is re- 
lated to some feature in the dark halo density profile, one could 
possibly expect to see some corresponding feature in the ro- 
tation curve or density profile at the radius at which the sur- 
face brightness profile starts to deviate from the fitted expo- 
nential disk. The radius which corresponds to the break points 
in the surface brightness profile have been indicated by verti- 
cal, dotted lines in the major-axis rotation curves (Fig.|5} and 
the density profiles derived from these (Fig. |8}. We do how- 
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Fig. 5. Radial /-band surface brightness profiles and B - V, V - I colour profiles derived for the target galaxies. The surface 
brightness profiles have been coiTected for inclination and Galactic extinction. The dashed lines represent the exponential disk 
profile determined from the outer isophotes. The radius at which the surface brightness profile starts to deviate substantially from 
an exponential disk is indicated by a vertical, dotted line. To allow a direct comparison between the linear and angular scales, 
arcseconds are used for the B - V profile and kiloparsecs for V - I. The radial axes of the graphs always extend to the same 
radius for each galaxy separately. Please note that the small-scale colour flucutations evident at large radii stem from decreasing 
signal-to-noise ratios in the outermost parts of the surface brightness profiles. To avoid cluttering, the error bars of the individual 
data points have not been included in this figure. 



ever not detect any obvious kinematic features associated with 
these breaks, except possibly for ESO 462-32, where the break 
coincides with the first data point of the plateau in the rotation 
curve. 

5. The central density parameter 

The theoretically predicted dark halo properties can be com- 
pared with observations in several different ways. One method 
is to examine the absolute value of the dark matter density in 
the inner regions of the dark halo. In this context, Alam et 
al. iZOOll suggested the use of a simple, dimensionless cen- 
tral density parameter Ay/i, defined as the mean dark matter 
density p (relative to the critical density of the Universe, pcrit) 
within the radius rv/2 where the galaxy rotation curve reaches 
half its maximum value Vmax: 

. Pjry/l) 1 / Vmax f 

Peril 2\Horv/2) 

Observationally, this quantity has the advantage that if the rota- 
tion curve rises at the outermost data point (which is common 
when only optical data is available), substituting V^ax for the 
outermost (highest velocity) point of the rotation curve will - 
for the range of likely density profiles - result in an upper limit 



on Av/2 (see Alam et al. 120021 for a discussion). Two other ef- 
fects also contribute in the same direction. By applying to 
the observed rotation curve without correcting for the rotational 
support provided by visible baryons (stars and gas), the central 
density of the dark halo will be overestimated. When compar- 
ing the observed and predicted dark halo Ay/i parameters, the 
observed value should also be lowered even further to correct 
for the halo contraction and density increase associated with 
baryonic cooling (e.g. Blumenthal et al. 119861 Gnedin et al. 
120041 Sellwood & McGaugh 2005i as current predictions for 
the dark halo Ay/i do not take this effect into account. 

Despite the fact that Ay/i will be overestimated when de- 
rived directly from observed rotation curves, previous inves- 
tigations have indicated that dwarf galaxies and LSBGs typ- 
ically display central halo mass densities significantly lower 
than those predicted for CDM halos in a ACDM cosmology 
(e.g. Alam et al. 2002 1. The scatter in Ay/i among the ob- 
served halos is also higher than expected from simulations (e.g. 
Hayashi et al. l2004at . 

In Fig. |5| we plot Vy^ax against Av/2 for LSBGs from de 
Blok et al. (PIMI, circles) and de Blok & Bosma dMHI trian- 
gles) together with our targets (filled squares). In all cases, the 
outermost data points of the average rotation curves have been 
taken as proxies for Vmax- The corresponding rv/2 have been 
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Fig. 6. The central density parameter Av/2 as a function of the 
maximum circular velocity. Markers indicate LSBGs from de 
Blok et al. (2001 ; circles), de Blok & Bosma (2002; triangles) 
and this paper (filled squares; ESO 146-14 in gray). The object 
F563-1, which appears in both the de Blok et al. and the de 
Blok & Bosma samples, has been indicated by filled markers. 
The lines represent the predictions of various models: standard 
ACDM (Alam et al. EUU^ thick soUd), ACDM with a tilted 
power spectrum and slightly different cosmological parame- 
ters (Alam et al. l2002l thin soUd), ACDM with additional 0.65 
eV neutrinos (Zentner & Bullock y2002i thick dashed), ACDM 
with a running mass index and n < 1 (Zentner & Bullock 2002 
thin dashed), CDM with a dark energy equation of state param- 
eter w = -1.5 (Kuhlen et al. 2005 , dash-dotted) and warm dark 
matter with a particle mass of 0.2 keV (Alam et al. 2002, dot- 
ted). Arrows indicate how our target galaxies would shift if a 
conservative coiTection for the baryonic disk was applied. For 
one object, ESO 548-09, this correction is too small to be seen. 
See main text for additional details. 

estimated by linear interpolation between the data points of the 
rotation curve. Due to its disturbed rotation curve, the position 
of ESO 146-14 (Vmax = 57 km/s and rv/2 = 3.1 kpc adopted) 
has been indicated by a gray marker. For this galaxy, the uncer- 
tainty mainly lies in the estimate of ry/i- Our Vmax estimate is 
consistent with that derived by Bergvall & Ronnback (1995 1, 
and using the width of the HI profile at 20% intensity instead 
of Ha data gives y^ax = 74 km/s (Mathewson & Ford 1996 1, 
which corresponds to only a minor shift of the data point un- 
less the rv/2 estimate is severely off. Interestingly, our galaxies 
are all located in the lower part of the observed Av/2 distribu- 
tion, where the conflict with theoretical predictions is the most 
severe. 

Included in Fig.|6lis the predicted Vmax-^v/i relation for a 
standard ACDM scenario (thick solid line; Alam et al. 2002) 
with a scale-free power spectrum of density fluctuations. Most 
of the data points are located at much lower halo densities than 
predicted by this model. The other lines represent various mod- 
ifications of standard ACDM suggested in the literature to im- 
prove the agreement with observations. Lower Ay/2 may be 
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Fig. 7. The dark halo central density parameter Av/2 as a func- 
tion of central surface brightness in the B-band. Markers in- 
dicate LSBGs from de Blok et al. dMHl circles), de Blok & 
Bosma ("2002" triangles) for which photometry is available. The 
object F563-1, which appears in both samples, has been indi- 
cated by filled markers. For our data, horizontal lines connect 
the true /Jsfi (filled squares; ESO 146-14 in gray) and that of an 
exponential disk extrapolated to the centre (open squares). See 
main text for additional details. 

achieved by tilting the power spectrum and adopting slightly 
non-standard values of the cosmological parameters Om, ^^a, 
erg and Ho (Alam et al. 2002 thin solid), adding a small con- 
tribution of hot dark matter in the form of 0.65 eV neutrinos 
(Zentner & Bullock .2002 thick dashed), by assuming a run- 
ning mass index at n < 1 (Zentner & Bullock 2002 thin dashed 
line) or by assuming that the dark energy has an equation of 
state (p - wp) characterized by w - -1.5 (Kuhlen et al. i20051 
dash-dotted line). The only model which is in reasonable agree- 
ment with our observations is one of the warm dark models 
presented by Alam et al. (120021 dotted line), where the dark 
matter is assumed to be a 0.2 keV fermion. It should be noted 
that the object ESO 120-021 from the sample of de Blok et al. 
(120019 . due to its very low Vmax, falls outside the range of the 
plotted model predictions in this diagram, and is therefore not 
included in the figure. 

In this comparison, we have adopted the minimum disk hy- 
pothesis, i.e. that the dark matter halo completely dominates 
the density of these galaxies at all radii. If this (unrealistic) ap- 
proximation is relaxed, our objects are expected to shift to even 
lower values of Ay/2. To demonstrate this, an estimated mini- 
mum correction for the contribution from the baryonic disk to 
the observed rotation curve has been indicated by arrows in 
Fig.|6l These corrections are derived from /-band data, assum- 
ing an exponential stellar disk with scale lengths given in Table 
121 In general, the existence of an empirical mass discrepancy- 
acceleration relation imposes strong constraints on the stellar 
population mass-to-light ratios of disc galaxies (e.g. McGaugh 
.2004 J . The colour-dependent mass-to-light ratios implied by 
this relation are, however, difficult to apply to our galaxies. 
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since both models (Portinari et al l2004> and empirical evidence 
(McGaugh l2005t indicates that the relation between colour and 
Ml Li, goes noisy and non-linear for objects as blue as these. 
Here, we instead adopt the smallest /-band stellar population 
mass-to-light ratio (M/L* - 0.7) found to be acceptable for 
the bluest LSBGs in the study of Zackrisson et al. ('2005X un- 
der the assumption of a Salpeter initial mass function. While 
HI maps are lacking for the galaxies in our sample, the con- 
tribution from neutral gas to the rotation curves of LSBGs is 
usually found to be smaller than that of the stellar disk within 
the optical rv/2. (e.g. de Blok et al. a 2001 , de Blok & Bosma 
12002 V We do however caution that the disk in many cases devi- 
ates substantially from the exponential form in the centre, and 
that this may affect the simple corrections applied here. 

In Fig. we plot Av/2 against the B-band central surface 
brightness levels yu^.o for those objects in Fig.|6lfor which pho- 
tometry is available. In the case of the de Blok et al. (2001 
circles) and de Blok & Bosma ('2002', triangles) data, the cen- 
tral surface brightness values correspond to that of an exponen- 
tial disk extrapolated to the centre. For these objects, the yu^.o 
values are in many cases based on /?-band data, converted un- 
der the assumption of B - R - 0.9 (following de Blok et al. 
1200 U . For our data, both the true, central yu^ o (filled squares) 
and that of the extrapolated exponential disk (empty squares) 
have been plotted. Despite a substantial scatter, the lack of 
faint ^Bfi objects at high Av/2 and bright iib.q objects at low 
Av/2 could possibly indicate the existence of a correlation be- 
tween these two parameters. While objects with widely differ- 
ent Vmax have been included in this plot, binning in Vmax does 
not seem to remove this correlation. There are at least two pos- 
sible interpretations of a relation of this kind. The first is that 
dark halos on average have much lower central densities than 
the ACDM scenario predicts, and that objects at extremely low 
central surface brightness levels simply allow a much cleaner 
measurement of this density, as the correction for baryonic ef- 
fects (which are usually assumed to increase the central den- 
sities) may be smaller for these objects. The other possibility 
is that extremely low surface brightness objects typically have 
dark halos with unusually low central densities, and that fo- 
cusing on such objects when comparing to the average results 
from ACDM simulations therefore introduces an unfair bias. 
These different possibilities are further discussed in Sect. 7. 
The choice of true or disk yu^ o in Fig. 0is not crucial for our 
objects, although our observations do not seem very remark- 
able in terms of surface brightness when disk fiBfl is plotted. 
In future investigations, it would be rewarding to explore how 
plots of central surface brightness versus Av/2 would look for 
objects with fixed Vmax and rotation curves corrected for the 
contribution from visible baryons. 

6. Slope of the central density profile 

Many attempts to test the predictions of the CDM model on 
the scales of galaxies have during recent years revolved around 
the so-called cusp/core problem. Whereas numerical simula- 
tions based on CDM have indicated that the central regions 
of dark halos should display a large increase in density - a 
cusp (e.g. NFW) - many observations have instead suggested 



the presence of a core of close to constant density (e.g. Blais- 



Ouellette et al. 120011 de Blo k & Bo sma 120021 Gentile et al. 
l2004l Spekkens & Giovaneni l2005i . The original NFW dai'k 
halo density profile is given by: 



p(/?) = 



Pi 



{RIR^){l+RIR5f' 



(2) 



where Rs represents a characteristic radius for the halo and 
Pi is related to the density of the Universe at the time of col- 
lapse. In the very centre, this density profile is characterized by 
p{r) oc R", where c w -1, i.e. a density cusp. The observations, 
on the other hand, often favour a value around a « 0, i.e. a 
constant-density core. From large samples of LSBGs and dwarf 
galaxies, both de Blok et al. (^003 1 and Spekkens & Giovanelli 
(12005 ) find on average a ^ -0.2. For this reason, cored density 
profiles have been advocated as superior dark halo models (e.g. 



Burkert ri995l de Blok et al. 2003 1. One commonly adopted 
such model is that of a pseudo-isothermal sphere, with a den- 
sity profile given by: 



piR) 



Po 



1 -1- 



(3) 



where po is the central density of the halo and Rq the radius of 
the core for which this density is representative. 

For several reasons, the slope of the central density pro- 
file may however not be the best test of CDM at the current 
time. On the theoretical side, there is still much controversy 
over the extent to which the NFW profile fitting formula re- 
ally is appropriate in the very centre. Whereas some authors 
have found central density slopes even steeper than the origi- 
nal NFW prediction (e.g. Fukushige & Makino 2003, Reed et 
al. 2005^1, others argue in favour of a more shallow slope (e.g. 
Ricotti Eb03l Stoehr 2006 1. On the observational side, a few 
studies have suggested that many measurements of the central 
density slope a may be biased by systematic effects so severe 
that observed cores may actually be consistent with intrinsic 
cusps (e.g. Swaters et al 2003; Rhee et al. '2004 ; Spekkens et 
al. 2005) - but see e.g. de Blok et al. (2003 1, de Blok (2004i 
and Gentile et al ( I2004> . for a different view. 

Despite these complications, it is nonetheless important to 
establish whether or not objects obeying our selection criteria 
deviate significantly from previous measurements of the cen- 
tral density profile. Under the assumption of a spherical matter 
distribution, the density profile can be calculated by direct in- 
version of the observed rotation curve V(Ry. 



1 LV(R)dViR) VHR) 



(4) 



This procedure assumes that the baryons in the disk act merely 
as test particles and have a negligible impact on the dynam- 
ics (the minimum disk hypothesis). The resulting density pro- 
files for the three objects (ESO 031-13, ESO 532-32 and ESO 
548-09) with the least irregular velocity fields in the innermost 
region are displayed in Fig.|8l 

To obtain the density slope a, a power-law is fitted to the 
innermost three data points of the density profile. Following de 
Blok et al 12003 ) . the adopted error bars on a correspond to the 
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Fig. 8. Density profiles derived under the assumption of a minimum disk and a spherical halo. Dashed lines indicate the power- 
law fitted to the innermost three data points of each profile. For comparison, the thick, dashed line indicates the forced fit of a 
cuspy density profile (a = -1) to the same data. The vertical dotted line marks the break radius at which the surface brightness 
profile starts to deviate from the slope derived from the outer isophotes. 



maximum change in slope introduced by including one data 
point more or one data point less in the fit. In Fig.|9l the result- 
ing central density slopes (large open markers) are compared 
to similar measurements (small black markers) by de Blok & 
Bosma (2002) and various models for the dark halo density 
profile (lines). Following current convention (e.g. de Blok & 
Bosma I2OO2I Spekkens & Giovanelh I2005> we plot a against 
Rin (representing the radius of the innermost point of the rota- 
tion curve), which gives an indication of the resolution of the 
observations. Since the rotation curve of ESO 548-09 becomes 
highly irregular at radii slightly outside that for which the den- 
sity slope was derived, this object has been indicated by a gray 
marker. Within the error bars, our measurements are in reason- 
able agreement with the de Blok & Bosma (20021 data and the 
average a ^ -0.2 found by de Blok et al. ( 2003 ) and Spekkens 
& Giovanelli (2005 1. They are furthermore in good agreement 
with the pseudo-isothermal sphere models (dashed lines, cor- 
responding to core radii of - 0.5, 1 and 2 kpc) advocated 
as models superior to the CDM predictions for the dark halos 
around galaxies. Due to the relatively large distances to the ob- 
jects in our samples, our observations are not particularly im- 
pressive in terms of resolution. They do however indicate cen- 
tral density slopes substantially shallower than those predicted 
by the original NFW profile (thin solid lines, corresponding to 
7?s = 3 and 30 kpc - suitable for halos in the dwarf-galaxy 
and large-galaxy mass range, respectively). Although the im- 
proved fitting formula (thick solid lines; same Rs) presented in 
Navarro et al. J2004t predicts a more shallow density slope in 
the very centre, the discrepancy at the radii probed by our ob- 
servations is not significantly smaller In fact, at the smallest 
radii where the simulations of Navarro et al. ( 2004 1 are con- 
sidered converged (vertical thin and thick dash-dotted lines for 
halos with R5 ~ 3 and 30 kpc, respectively) the agreement be- 
tween the predictions (a x -1.4 and -1.2, respectively) and 
the observations is worse than when the original NFW profile 
is used. 



7. Discussion 

The very low central mass concentrations (as measured by 
Av/2) of the halos around these blue LSBGs appear to be in 
serious conflict with the predictions of the halos formed in 
the ACDM scenario. We furthermore find no sign of the steep 
(a < -1) CDM central density cusps, which are still advocated 
by some authors (e.g. Fukushige & Makino 2003^ Navarro et 
al. ETO Reed et al. l2005l . 

Since those of our galaxies for which the inner density pro- 
file has been estimated appear to be completely bulgeless, the 
lack of the high central concentrations predicted for CDM ha- 
los contradicts the mechanism suggested by Mo & Mao J2004> 
for lowering the central CDM halo density, as this model sug- 
gests that the original CDM halo profile should still be valid 
for bulgeless galaxies. Although only optical surface bright- 
ness profiles are presented in this paper, Bergvall et al. ( I1999> 
have shown that the decrease in surface brightness towards the 
centre (compared to an exponential disk profile), which is here 
interpreted as the absence of a bulge, typically persists in the 
near-lR as well. Hence, it cannot easily be attributed to dust 
effects. 

Before claiming a general failure of the CDM paradigm, a 
number of additional complications do however require con- 
sideration. 

7.1. Issues related to disk inclination 

The inclinations of the galaxies in our sample have been es- 
timated to lie in the range / - 76-83°. Although selecting 
high inclination disks allows the study of galaxies at fainter 
surface brightness levels than would otherwise have been pos- 
sible, this particular property of the galaxies studied also opens 
the possibility that conclusions drawn from the measured ro- 
tation curves may be affected by projection effects. Matthews 
& Wood ( .2001 J modelled the effects of internal extinction and 
projection effects in the disks of LSBGs, concluding that ro- 
tation curves should be negligibly affected at inclinations of 
/ < 85° even at high optical depths. Some effects could however 
be discerned even at / - 85° in the case of a very steep rota- 
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tion curve slope. Similar calculations were carried out by Baes 
et al. J2003> . showing that projection and dust effects are ex- 
pected to be very small, although not undetectable, at / < 85°. 
The optimal inclination for accurate rotation curve retrieval ap- 
pears to be around / w 60°. S waters et al. (2003 1, on the other 
hand, simulated the various inclination-related systematic ef- 
fects involved in measuring the central density profile, with the 
alarming result that measurements of the inner slope of the den- 
sity profile could be significantly biased even in galaxies with 
; ~ 60°. These results were later questioned by de Blok et al. 
J2003> . who found the observed inner slope distribution to be 
insensitive to whether / < 80° or ; < 85° was adopted in the 
galaxy samples used. Rhee et al. J2004> nonetheless detected 
an inclination-related bias of the type envisioned by Swaters et 
al. in the de Blok et al. data. 

Since the line of sight traverses the most diverse orbits of 
the observable disk at small distances from the centre, projec- 
tion effects are however expected to become smaller further out 
(e.g. Rhee et al. 120041. The concentration parameter c (often 
defined as the ratio between the virial radius and the character- 
istic halo radius of the NFW profile: c = RyiJRs) in many ways 
provides a more robust test than the innermost slope of the den- 
sity profile, since it depends on the rotation curve behaviour 
over a much larger radial interval. Figs. 8 & 9 of Swaters et al. 
J200 3 1 do indeed confirm that the inclination-dependent sys- 
tematic effects should be smaller for c than for the innermost 
slope of the density profile. The related central density param- 
eter Av/2 used here also extracts information from much larger 
radii, and should be less affected by inclination-related effects. 

In summary, a certain degree of scepticism should be ap- 
plied when interpreting the inner slope of the density profile 
measured in high-inclination galaxies. The patchy distribution 
of Ha emission seen at the very faint surface brightness levels 
of the galaxies analyzed here could further augment this prob- 
lem. To put the conclusions reached in this paper on a more 
robust footing, it would therefore be very useful to extend this 
study to a sample of galaxies selected with similar colour and 
surface brightness criteria but lower inclinations. 

7.2. Issues related to the use of long-slit spectroscopy 

Aside from the extinction and projection effects that complicate 
the interpretation of rotation curves of high-inclination disk 
galaxies, there are also other observational effects that could 
affect the dark halo density profiles inferred from long-slit ob- 
servations. Wide slits, bad seeing and slits not properly posi- 
tioned along the major axis would all make the intrinsically 
cuspy density profiles appear more core-like, as demonstrated 
by the simulations of de Blok et al. (^2003 1 and Swaters et al. 
J2003> . With good seeing (here < 0.8") and narrow slits (1"), 
the first two effects are however predicted to be modest. While 
slit offsets can certainly be a problem for individual objects 
with complicated isophotes (e.g. ESO 546-34 in our sample), 
the typical position eiTors of our observations are estimated to 
be no more than 1", whereas offset errors of 3^" would be re- 
quired for intrinsically cuspy halos to mimic cores. 
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Fig. 9. Central slope of the dark halo density profile (a) vs. 
the radius of the innermost point in the rotation curve, Ri^. 
Markers indicate the measurements of a from our sample (large 
white/gray markers) and measurements from the LSBG sam- 
ple of de Blok & Bosma (2002 small black markers). Dashed 
lines indicate the density slopes of pseudo-isothermal halos 
with Rc = 0.5, 1 and 2 kpc (from left to right). Thin solid 
lines represent the density slopes of NFW halo profiles with 
/?s = 3 and 30 kpc (from bottom to top in the rightmost part of 
the plot). Thick solid lines represent the corresponding slopes 
from the updated Navarro et al. 12004) halo profile. The verti- 
cal thin (/?s ~ 3 kpc) and thick (Rs « 30 kpc) dash-dotted lines 
indicate the innermost radii at which the Navarro et al. (I2004> 
simulations are considered converged. 



7.3. Issues related to the shapes of dark matter halos 

The analysis presented in this paper is based on the assumption 
of a spherical dark halo, while the CDM scenario in fact pre- 
dicts dark halos to be triaxial. In Hayashi et al. ( 2004b i, it was 
argued that the discrepancy between observed LSBG rotation 
curves and the predicted properties of CDM halos could be re- 
duced by properly taking the complicated dynamical effects of 
a disk rotating in a triaxial dark halo potential into account. As 
demonstrated in Sect. 3, our minor-axis data do however not 
show any sign of the kinematical signature derived by them for 
the elliptical disk that formed in their simulation. This is con- 
sistent with the results from other studies which also suggest 
that the observed non-circular motions in disks are insufficient 
to mask central density cusps (Gentile et al. 2005 1 and that the 
elUpticities of galactic disks in general are small (see Combes 
120021 and references therein). Even so, this does not rule out 
the possibility that the apparent discrepancy between our re- 
sults and the CDM predictions may become smaller once the 
approximation of a spherical halo is relaxed. As Hayashi et al. 
( I2004b> present only a single viewing angle from a single sim- 
ulation, for which a particular orientation of the disk inside the 
triaxial dark halo was assumed, it is not at all clear how generic 
their results really are. 
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Recently, Dutton et al. J2005> investigated the effect of 
relaxing the assumption of spherical halos in favour of axi- 
ally symmetric ellipsoids when interpreting rotation curve data, 
and found the best-fit halo concentration parameter to decrease 
for a halo flattened in the direction perpendicular to the disk, 
and to increase for a halo elongated in the same direction. As 
previous attempts to fit NFW halos to LSBG rotation curve 
data have indicated that the average concentration parameter 
c tends to be too low compared to the CDM predictions (e.g. 
McGaugh et al 2003 de Blok et al. 2003 1, a flattened halo 
would hence increase the tension between CDM predictions 
and LSBG observations, whereas elongated halos would de- 
crease it. Dissipationless CDM simulations (e.g. Jing & Suto 
I2002> predict most halos to be prolate, i.e. having two axes of 
similar size and a third axis that is longer - which is inciden- 
tally the halo shape adopted in the Hayashi et al. simulation 
as well. However, the typical halo shape can be altered sub- 
stantially by baryonic processes (e.g. Kazantzidis et al. 2004) 
and may even become oblate (e.g. Dubinski 1994j, i.e. having 
two axes of similar size and third axis which is shorter. Indeed, 
many observational investigations (see Sackett '1999' for a re- 
view) suggest that the dark halos of galaxies truly are oblate. As 
LSBGs are often assumed to be associated with dark halos of 
unusually high spin (e.g. Boissier et al. '2003V and simulations 
indicate a correlation between high spin and oblateness (Moore 
et al. 2004 1, it seems quite reasonable to assume that LSBGs 
are located inside oblate halos. As hydrodynamic simulations 
of disks forming in dark matter halos furthermore indicate that 
the rotational axis of the disk tends to align with the minor 
axis of the inner halo (Bailin et al. I2005bl >. it also seems rea- 
sonable to assume that an oblate halo should coiTespond to the 
polar flattening scenario investigated by Dutton et al. (2005 1. 
This suggests that relaxing the assumption of a spherical halo 
would just strengthen the case against CDM even further. There 
is, however, at least one potential flaw in this argument. While 
simulations predict a relation between the angular momentum 
of a dark halo and its shape, they also predict a relation between 
angular momentum and the concentration parameter c (Bailin 
et al. 1200 5a I. which could bias measurements of c in LSBGs 
towards values lower than the cosmic average (as further dis- 
cussed in Sect. 7.5). 

To summarize, it is very difficult to say anything conclu- 
sively about how relaxing the assumption of a spherical dark 
halo would affect the comparison between CDM predictions 
and observations of LSBGs. To settle these issues, both im- 
proved predictions for the shapes of CDM halos under the in- 
fluence of baryonic processes and methods for analyzing the 
kinematic data in the framework of these non-spherical halos 
would be required. 

7.4. Issues related to the possible influence ofbaryons 

The reason for using LSBGs to test CDM halo predictions is 
that the effects of baryons on the overall density profile are be- 
lieved to be minimized in these systems. Because of the strong 
correlation between dynamical mass-to-light ratio of and cen- 
tral disk surface brightness, the faintest LSBGs would appear 



to be the best targets available. Although the presence of un- 
detected disk components in LSBGs have been claimed (Fuchs 
I2003I I. this does not challenge the fact that in LSBGs, the ra- 
dius at which the dark matter halo dominates over the maxi- 
mum disk by a factor of two is quite small (McGaugh & de 
Blok'1998aV although with a substantial scatter. More wor- 
risome is the claim of Graham (2002), that the McGaugh & 
de Blok relation may arise from selection effects, and that all 
LSBGs need not be as dark-matter dominated as previously as- 
sumed. Indeed, Bizyaev & Kajsin (2004) find that for LSBGs 
with bulges, the ratio of dark to luminous matter is not sub- 
stantially different from that of high surface brightness galax- 
ies. If the inner regions of LSBGs are not as CDM dominated 
as previously believed, but instead dominated by dark baryons 
(e.g. Combes 2004), this opens the possibility that some bary- 
onic process - e.g. bar formation (Athanassoula 2004) - may 
have altered the inner density profile in some complicated way 
not taken into account in the simple scenario of adiabatic con- 
traction (Blumenthal et al. 1986 Gnedin et al. 2004; Sellwood 
& McGaugh 2005), which is often advocated to argue that 
baryons should increase - and not decrease - the central den- 
sity. The commonly adopted procedures for analysing velocity 
profiles of disk galaxies can also lead to underestimates of both 
the central density and the slope of the central density profile 
if undetected bars are present (Rhee et al. .2004j . These issues 
need to be investigated further. 

7.5. Issues related to dark halo selection effects 

Finally, there is the issue of whether or not the dwarf galax- 
ies and LSBGs which are currently used to test CDM pre- 
dictions on galactic scales constitute a biased dark halo sam- 
ple. If LSBGs are preferentially formed in the low-density tail 
of the halo distribution, then discrepancies of the type high- 
lighted here may perhaps not pose any serious threat to the 
CDM scenario (e.g. Zentner & Bullock [20021 Jimenez et al. 
2003 Bailin et al. 2005a), since comparing the dark halo prop- 
erties of the galaxies with the lowest surface brightness lev- 
els to predictions for typical CDM halos would then be quite 
unfair. A bias of this kind can for instance arise if the halos 
surrounding LSBGs preferentially have an unusually late for- 
mation time (e.g. Jimenez et al. I2003t or if LSBGs are asso- 
ciated with dark halos of unusually high spin (e.g. Boissier et 
al. l2003> . Due to the degeneracy between age and star forma- 
tion history, proving that LSBGs are young systems is how- 
ever not trivial (Zackrisson et al. 2005 ). It is furthermore not 
obvious that the cosmic distribution of halo spin is consistent 
with the observed surface brightness distribution of galaxies. 
In the model of Boissier et al. (I2003t . only a minority of disk 
galaxies (^ 35%) become LSBGs, which does not appear to be 
consistent with current estimates of the size of the LSBG pop- 
ulation. With the LSBG definition adopted by Boissier et al. 
(jUo.B Z 22 mag arcsec"^), the surface brightness distribution of 
disk galaxies presented by O'Neil & Bothun ( I2000> would for 
instance imply that > 80 % are LSBGs. McGaugh & de Blok 
tl998 a) furthermore argue that the notion that LSBG reside in 
halos of unusually low density would be difficult to reconcile 
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with the observed Tully-Fisher relation. Clearly, more detailed 
studies of halo concentration bias as a way to reconcile LSBG 
observations with CDM halo predictions are urgently needed. 

8. Summary 

By studying the Ha kinematics of a sample of LSBGs with 
unusually blue colours and very low surface brightness centres, 
we find that: 

- These galaxies all fall inside the region of the diagram of 
maximum velocity V^ix vs. central density parameter Ay/i, 
which poses the toughest challenge for current ACDM dark 
halo predictions. The predictions of a number of alternative 
dark matter and dark energy models are also shown to be in 
conflict with these data. 

- By combining our observations with data from the litera- 
ture, we uncover a possible relation between central surface 
brightness jUb o and the halo central density parameter Av/i- 
We argue that such a relation may explain the previous re- 
sult and suggest that the faintest LSBGs - in general - may 
provide the most discriminating tests of ACDM on galactic 
scales, unless these galaxies systematically tend to form in 
the low-density tail of the dark halo distribution. 

- Our estimates of the slope of the innermost density profile 
in the centre of these galaxies are more reminiscent of con- 
stant density cores (consistent with previous investigations 
of other LSBGs and dwarf galaxies) than the steep density 
slopes predicted for CDM halos by either the original NFW 
fitting formula or the improved formula by NavaiTO et al. 
(120041 . 

- Since most of the galaxies analysed are bulgeless, our in- 
ability to confirm the CDM halo predictions disfavours the 
scenario suggested by Mo & Mao ( I2004> . in which the halo 
density profile is altered by feedback associated with a fast 
collapse phase, as this model predicts that dark halos with 
the original CDM halo profile should still be detectable in 
bulgeless galaxies. 

- Our minor-axis data show no resemblance to the kinemati- 
cal signature found by Hayashi et al. J2004b> in their sim- 
ulation of a disk rotating in a triaxial CDM halo. 

- The inner radius at which the surface brightness profile 
starts to deviate from an exponential disk derived from the 
outer isophotes does not seem to correlate with any feature 
in the rotation curve or the estimated, spherically averaged 
density profile. The origin of these features remains a mys- 
tery. 
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